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We analyze the theory of optimal investment in risky assets, developed recently by Marsili, Maslov 
and Zhang [Physica A 253 (1998) 403]. When the real data are used instead of abstract stochastic 
process, it appears that a non-trivial investment strategy is rarely possible. We show that non-zero 
transaction costs make the applicability of the method even more difficult. We generalize the method 
in order to take into account possible correlations in the asset price. 
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I. INTRODUCTION 

Non-equilibrium statistical mechanics, especially the 
theory of stochastic processes, finds recently wide appli- 
cability in economics. First area, intensively studied in 
the last several years, is the phenomenology of the signal 
(price, production, and other economic variables) mea- 
sured on the economics system Scaling concepts 
proved to be a very useful tool for such analysis. 

Second area concerns optimization. In the competitive 
economics, agents should maximize their survival proba- 
bility by balancing several requirements, often mutually 
exclusive, like profit and risk |p|-|l3||. Third area com- 
prises creation of models which should grasp particular 
features of the behavior of real economics, like price fluc- 
tuations [||-^9). 

We focus here on an aspect of optimization, discussed 
recently by Marsili, Maslov and Zhang pfjfl . In a sim- 
plified version of the economy, there are two possibilities 
where to put a cash: to buy either a risky asset (we shall 
call it stock, but it can be any kind of asset) or a risk- 
less asset (deposit in a bank). In the latter case we are 
sure to gain each year a fixed amount, according to the 
interest rate. On the contrary, putting the money en- 
tirely to the stock is risky, but the gain may be larger 
(sometimes quite substantially). We may imagine, that 
increasing our degrees of freedom by putting a specified 
portion of our capital into the stock and the rest to the 
bank may lead to increased growth of our wealth. This 
way was first studied by Kelly and followers |^l|,^2| and 
intensively re- investigated recently 1 20 2^-|27[| . 

The point of the Kelly's approach is, that if we sup- 
pose that the stock price performs a multiplicative pro- 
cess [^8-31 , the quantity to maximize is not the average 
value of the capital, but the typical value, which may 
be substantially different, if the process is dominated by 
rare big events. It was found that given the probability 
distribution of the stock price changes, there is a unique 
optimal value of the fraction of the investor's capital put 
into the stock. 

The purpose of the present work is to investigate the 



practical applicability of the strategy suggested in |20|,|23| . 
Let us first briefly summarize this approach. We suppose 
that the price pt of the stock changes from time t to t + 1 
according to a simple multiplicative process 



Pt+i = Pt e Vt 



(1) 



where r\t for different t are independent and equally dis- 
tributed random variables. The angle brackets <> will 
denote average over these variables. 

We denote Wt the total capital of the investor at the 
moment t. The fraction r of the capital is stored in stock 
and the rest is deposited in a bank. We will call the 
number r investment ratio. The interest rate provided 
by the bank is supposed to be fixed and equal to p per 
one time unit. The strategy of the investor consists in 
keeping the investment ratio constant. It means, that 
he/she sells certain amount of stock every time the stock 
price rose and sell when the price went down. 

If we suppose that the investor updates its portfolio 
(i. e. buys or sells some stock in order to keep the in- 
vestment ratio constant) at each time step, then starting 
from the capital Wq, after N time steps the investor owns 
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(2) 



The formula can be simply generalized to the situation 
when there is a non-zero transaction cost equal to 7 (see 
also |27)) and the update of the portfolio is done each M 
time steps. We assume for simplicity that TV is a multiple 
of M. 
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(1 + p) M + r(e^» (1 + G) - (1 + p) M ) 
1 + rG 



(3) 



where we denoted r\Mt = X)f=M* f 1 7 li an d G — 
7 sign(M ln(l + p) - fj M t) ■ 

We can see that like the stock price itself, the capital 
performs a multiplicative process. 
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et(r)W t 
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where the random variables et(r) depend on the invest- 
ment ratio as a parameter. 

For N sufficiently large the typical growth of the cap- 
ital (Wt+i/Wf) typical is not equal to the mean < e(r) > 
as one would naively expect, but is given by the median 
pQ|, which in this case gives 



7(r) = \og{{W t+ i/Wt)ty P icai) =< loge(r) > 



(5) 



Therefore we look for the maximum of g as a function 
of r, which in the simplest case without transaction costs 
leads to the equation 
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l + p + r opt (c f '-l-p) 



>= 



(6) 



for the optimum strategy r opt . If the solution falls outside 
the interval [0, 1], one of the boundary points is the true 
optimum, based on the following conditions. If g'(0) < 
the optimum is r op t = 0. If g'(l) > the optimum is 
r op t = 1. 

If r\ is a random variable with probability density 
1 



(5(rj — m — d) + 5(j] — m + d)) 



(7) 



moment, but not necessarily). Then, we imagine that 
at time t\ an investor started with capital Wt 1 — 1 and 
during that period followed the strategy determined by 
certain value of r. We compute his/her capital W t2 (r) 
at final time and find the maximum of the final capital 
Wt 2 (r) with respect to r. We call the value r opt maximiz- 
ing the final capital two-time optimal strategy. Optimum 
strategy in the past can be then used as predicted opti- 
mal strategy for the future. 
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FIG. 1. Time evolution of the NYSE composite index. 
Time is measured in working days from t = which is 2 
January 1990 to t = 2181, which is 31 December 1998. The 
vertical axis is in logarithmic scale. 



the solution of (|6J) is straightforward: 
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(8) 



In more complicated cases we need to solve the equa- 
tion ^ numerically. However, for small mean and vari- 
ance of 1] approximative analytical formulae are fairly 
accurate [ p3|p^ ]. We found, that equally good approx- 
imation is obtained, if we set m —< rj > and d = 
\f< ri 2 > - < rj > 2 in the Eq. |. 

In the next sect ion we investigate the method with real 
data. 



Section [III 
costs. In Sec 



shows the influence of the transaction 
1V| a generalization of the method for the 
case of time-correlated price is shown. Finally, in Sec. |v| 
we discuss the obtained results. 



II. TWO-TIME OPTIMAL STRATEGIES 

In the previous section we supposed the following pro- 
cedure: the investor takes the stock price data and ex- 
tracts some statistical information from them. This infor- 
mation is then plugged into theoretical machinery, which 
returns the suggested number r. However, we may also 
follow different path, which should be in principle equiv- 
alent, but in practice it looks different. 

Namely, suppose we observe the past evolution of the 
stock price during some period starting at time t\ and 
finishing at time ti (most probably it will be the present 



The capital at time ti is again 



W t2 (r) = Y[(l + p + r(ev -1-p)) 



(9) 



t=ti 



and its maximization with respect to r leads to equation 
t 2 -i 



g'{r opt ) = ^2 i 



e'" 



t=ti 



+ p + r opt (c"* - l-p) 



(10) 



which gives the optimal investment ratio J" op t(ti, ta) as a 
function of initial time t\ and final time ti. Note that it 
is an analog of the equation (g) but we deal with time 
averages here, not with sample averages as before. This is 
also another justification of the procedure of maximizing 
< \og(W t +i/W t ) > instead of < W t +i/W t >■ 

For comparison with reality we took the daily values 
of the New York Stock Exchange (NYSE) composite in- 
dex. The time is measured in working days. The period 
studied started on 2 January 1990 (t = 0) and finished 
on 31 December 1998 (t = 2181). The time evolution of 
the index x(t) is shown in Fig. |l|. The values of rj are 
determined by exp(7y t ) = x(t + l)/x(t). 

The data of NYSE composite index were analyzed by 
calculating the two-time optimal strategies r op t(ii,£2)- 
As a typical example of the behavior observed, for ini- 
tial time ti — 300 we vary the final time £2 up to 2180. 
We used the interest rate 6.5% per 250 days (a realistic 
value for approximately 1 year). In this case we neglect 
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the transaction costs, 7 = 0. The influence of non-zero 



transaction costs will be investigated in Sec. 111. The re- 
sults are in Fig. ||(a). We investigated also the possibility 
that the investment ratio goes beyond the limits and 1 , 
which means that the investor borrows either money or 
stock. We imposed the interest rate 8% on the loans and 
calculated again the optimal r. The results are in Fig. 
||(c). We can see several far-reaching excursions above 
1 and some also below 0, which indicates that quite of- 
ten the optimal strategy requires borrowing considerable 
amount of money or stock. 
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FIG. 2. The two-time optimal investment ratio for interest 
rate 6.5% per 250 days. The initial time is 300. The trans- 
action costs are 7 = (a) and 7 = 0.005 (b). In (c), loans 
are allowed with interest rate 8% per 250 days, transactions 
costs are 7 = 0. 

An important conclusion may be drawn from the re- 
sults obtained: the optimal strategy r op t(£i, t%) as a func- 
tion of the final time ti does not follow any smooth tra- 
jectory. On the contrary, the dependence is extremely 
noisy, as can be seen very well in the Fig. ||. Moreover, 
the strategy is very sensitive to initial conditions. If we 
compare the strategy r op t(ti, £2) and r op t(ii + Ai, £2) for 
slightly different initial time, big differences are found in 
regions, where the strategy is non-trivial (0 < r opt < 1). 
In Fig. H we show for At = 1 the average difference in 
optimal strategy 

Ar opt (t) = (|r op t(ti,*i+t)-r <v t(ti + l,*i + t)|) (11) 

where the average is taken over all initial times t\ with 
the constraint, that we take into account only the points 
where both optimal strategies r op t(ii, t\ +t) and r op t(ii + 
1, ti + t) are non-trivial. 



FIG. 3. Average difference in optimal strategy when the 
initial times differ by 1 day. Only points where the strategies 
are non-trivial are taken into account. 

Due to poor statistics, the data are not very smooth. 
We can also observe apparent two branches of the depen- 
dece, which is caused by superimposing data from differ- 
ent portions of the time evolution of the index. However, 
despite of the poor quality of the data, we can conclude, 
that even after a period as long as 1000 days (approx- 
imately 4 years) the difference of 1 day in the starting 
time leads to difference in optimal strategy as large as 
about 0.2. This finding challenges the reliability of the 
investment strategy based on finding optimal investment 
ratio r. 

Moreover, we can see that if loans are prohibited, there 
are long periods where the optimal strategy is trivial 
(''opt = or r opt = 1). We investigated the whole his- 
tory of the NYSE composite index shown in Fig. [I] and 
determined, for which pairs (£1,^2) the optimal strategy 
r pt(ti,t2) is non-trivial. In the Fig. || each dot repre- 
sents such pair. (In fact, not every point was checked: 
the grid 5x5 was used, i. e. only such times which are 
multiples of 5 were investigated.) 

We can observe large empty regions, which indicate ab- 
sence of a non-trivial investment. In order to understand 
the origin of such empty spaces, let us consider a simple 
model. Suppose we have the random variable distributed 
according to (^), and p = 0. Then the conditions for the 
existence of non-trivial optimal strategy between t\ = 
and t2 = N are 
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?'(0) = £ (e* - 1) > 



and 
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7 '(l) = X;(l-e-*)<0 



(12) 
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FIG. 4. Existence of non-trivial strategies: each dot rep- 
resents starting and final time between which a non-trivial 
optimal strategy is found. 

Let us compute the probability p nt that both of these 
conditions are satisfied. We have 



g'(Q) = N{e m coshd- 1) +c f 
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' sinh d z t 



t=o 



and 



g'(l)=N(l 



cosh d) 



N-l 

' sinh d Zt 



(14) 



(15) 



t=o 



where z's can have values +1 or -1 with probability 1/2. 



The sum X^cj 1 z * nas binomial distribution, and for 
large N we can write 

"VjV (coth d—c m / sinh 



Pnt 



dC 



//V(cothd-G- m /sinh d) 



exp(-£) . (16) 



We can see immediately that p n t has a value close to 
1 for the number of time steps at least 



N 



(17) 



For the data in Fig. [l] we found d ~ 0.01, which means 
N ~ 10000 days, or 40 years. This is thus an estimate 
of how long we need to observe the stock price before 
a reliable strategy can be fixed. However, during such a 
long period the market changes substantially many times. 
That is why no simple strategy of the kind investigated 
here can lead to sure profit. 



III. TRANSACTION COSTS 

We investigated the influence of the transaction costs 
7 and time lag M between transactions. We found nearly 
no dependence on M, but the dependence on 7 is rather 
strong. It can be qualitatively seen in Fig. ||(b). When 
we compare the optimal strategy for 7 = and 7 = 0.005, 
we can see, that already transaction costs 0.5% decrease 



substantially the fraction of time when the strategy is 
non-trivial. We investigated the dependence of the frac- 
tion /nontriviai of time pairs {ti,t 2 ) between which a non- 
trivial strategy exists on the transaction costs. We have 
found that it decreases with 7 and reaches negligible 
value for 7 > 0.006. This behavior is shown in Fig. || 
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FIG. 5. The dependence of the fraction of time pairs, be- 
tween which a non-trivial investment optimal strategy exists, 
on the transaction cost. The time interval investigated is from 
time to time 1600. 

The explanation of this behavior lies in the fact, that 
the transaction costs introduce some friction in the mar- 
ket, which means that large changes of the investment 
ratio are suppressed. Because the investment ration is 
mostly or 1 even for 7 = 0, this implies that changing r 
from or 1 to a non-trivial value is even harder for 7 > 
and a non-trivial investment becomes nearly impossible 
for large transaction costs. 



IV. INVESTMENT IN PRESENCE OF 
CORRELATIONS 

In order to improve the strategy based only on the 
knowledge of the distribution of 77, we would like to in- 
vestigate a possible profit taken from the short-time cor- 
relations. 

Imagine again the simplest case, when 77 can have only 



two values, r\ 



+ — 



-d and 77 = m — d. However, now r/ t 



and r\t-\ may be correlated and we suppose the following 
probability distribution P(j] t -i, rj t ) — 1/4 + c if rjt-i = rjt 
and P(r)t-i,T)t) — 1/4 — c if 774-1 7^ i] t . The parameter c 
quantifies the short- time correlations. 

At time t the strategy r(r] t _ 1 ) should depend on the 
value of 77 in the previous step. In our simplified situation 
we have only two possibilities, r + = r(rj + ) and r~ = 
r(?7 _ ). The problem then reduces to maximization of the 
typical gain 

g(r+, r~) =< ln(l + p + r(77 t _ 1 )(e" t - 1 - p)) > (18) 



which leads to decoupled equations for r^ pt and 
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(19) 
(20) 

(21) 
(22) 



The solution is a straightforward generalization of Eq. 

The above procedure works equally well even in the 
case of more complicated time correlations. For example 
we may imagine that the price evolution is positively cor- 
related over two time steps, i. e. Prob(?7t_2 = r)t) > 1/2, 
while Prob(?7i_i = r) t ) — 1/2. Generally, we have some 
joint probability distribution for the past and present 
P(r/ < ,?y), where we denote r/< = [..., r]t-3, r) t -2, Tjt-l] and 
r\ = r\ t . The typical gain becomes a functional depending 
on the strategy r(?7 < ) which itself depends on the past 
price history. 

However, maximizing this functional by looking for its 
stationary point leads to very simple set of decoupled 
equations for the strategies 



d?7 P(r? < ,r?) 



-l-p 



l + p + r opt (rj<)(e"-l-p) 



= . (23) 



In the simplest case, when we assume that the strategy 
depends only on the sign of r\ in the previous step, we per- 
formed the analysis on the NYSE composite index shown 
in the Fig. |]. We found optimal pairs [r„ t , r~ t ]. Con- 
trary to the case when correlations were not taken into 
account, no non-trivial investment strategy was found. 
So, instead to improve the method of Ref. p0|] , the ap- 
plicability of this method is further discredited. 



V. CONCLUSIONS 

We investigated the method of finding the optimal 
investment strategy based on the Kelly criterion. We 
checked the method on real data based on the time evo- 
lution of the New York Stock Exchange composite index. 
We found, that it is rarely possible to find an optimal 
strategy which would be stable at least for a short pe- 
riod of time. There are several reasons, which discredit 
the method based on the Kelly criterion. First, the op- 
timal investment ratio fluctuates very rapidly in time. 
Second, it depends strongly on the time, when the in- 
vestment strategy started to be applied. The difference 
of 1 day in the starting moment makes substantial differ- 
ence even after 1000 days of investment. Third, the frac- 
tion of days, for which a non-trivial investment strategy 
is possible, is very low. This fraction also decreases with 



transaction costs and reaches negligible values for trans- 
action costs about 0.6%. Taking into account possible 
correlations in the time evolution of the index makes the 
situation even less favorable, reducing further the frac- 
tion of times, when a non-trivial investment is possible. 

We conclude, that straightforward application of the 
investment strategy based on the Kelly criterion would 
be very difficult in real conditions. The question re- 
mains, whether there are other optimization schemes, 
which would lead to more certain investment strategies. 
It would be also interesting to apply the approach used in 
this paper in order to check the reliability of the option- 
pricing strategies. 
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